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We investigate the statistical properties of a randomly branched 3-functional TV-link polymer 
chain without excluded volume, whose one point is fixed at the distance d from the impenetrable 
surface in a 3-dimensional space. Exactly solving the Dyson-type equation for the partition function 
Z{N,d) = N- 9 e lN in 3D, we find the "surface" critical exponent $ = ~, as well as the density 
profiles of 3-functional units and of dead ends. Our approach enables to compute also the pairwise 
correlation function of a randomly branched polymer in a 3D semi-space. 



I. INTRODUCTION 



The rapid development of statistical physics of macromolecules is indebted, in first turn, to the representation of 
linear polymers by Markov chains. The application of a theory of Markov processes for a description of conformational 
properties of linear polymer chains has become invaluable Besides the linear chains, the polymers of complex 

geometry, such as, for example, stars, brushes and branched chains, are of extreme importance in physical and chemical 
properties of macromolecular compounds 0, 0, 0] ■ 

Randomly branched macromolecules are ones of the most interesting polymer systems with nonlinear structure. The 
possibility to form topologically different architectures leads to an extra contribution to the conformational entropy 
of samples constituted by randomly branched chains. There is a deep relation between physics of randoml y b ranched 
pol ym ers and "lattice animals" appeared in percolation 0,0,01 and that of the gelation processes [fj HoL TTTL Il2l IT^L 
H3.lla.llq . There are various methods for statistical description of annealed randomly branched chains. On the one 
hand, we know successful attempts to compute various conformational properties of randomly branched polymers in 
terms of branched Markov processes 0, 0, 0, . 

Similarly to the theory of linear polymer chains (see, for example, 0, 0|), the problem of taking the volume 
interactions (i.e. interaction of particles situated topolgically far from each other) into account plays the central role 
in the construction of proper theory for the randomly branched chains. Concerning this problem, one can imagine 
two most interesting limiting cases, corresponding to strong repulsive volume interactions and negligibly small ones, 
respectively (we discuss the applicability of these two limiting cases in the last section of this paper). 

Speaking of the limit of strong repulsive volume interactions, the special attention should be paid to the application 
of the supersymmetry in the description of thermodynamic characteristics of randomly branched chains. During last 
two decades since the seminal works [2(| HU (see also the short self-reviews |2^|) the sypersymmetry becomes 
a powerful tool for the investigation of statistics of randomly branched chains. In particular one has to note the 
recent contributions |3Ct l3l| which give the impact in understanding the very origin of the supersymetric dimensional 
reduction for randomly branched polymers. Speaking of the main objective of the present work - the influence of 
boundary conditions on statistical properties of randomly branched polymers, we know some recent successful attempts 
of calculating the critical exponents of the randomly branched chains with volume interactions near the impenetrable 
wall and the thorough study of absorption of a randomly branched polymer |3^| . The former of these works is 
based on rather general lattice approximation (note that lattice theories always correspond to the presence of excluded 
volume), the latter is based on the supersimmetry technique. In the discussion section we will compare the results of 
32, 331 with ours obtained in the limit of no volume interactions. 



In the no interaction (ideal) case the basic question of the theory concerns the evaluation of the partition function 
Z(N) of a non-self-interacting branched structure without loops. Briefly, the problem is formulated as follows. Take N 
elementary units (monomers) such that each monomer has no more than f max branches (we call them / maa; -functional 
monomers) and compute all possible ways to arrange them in a single- connected loopless cluster. We call such a cluster 
"a / maa; -branching random tree" . In what follows we shall consider for simplicity only the case f m ax = 3. The problem 
of calculating the partition function Z(N) in the free space does not meet any difficulties and has been solved by 
various methods (see, for example, |22u23|). However the influence of boundary conditions on statistical properties of 
randomly branched polymers is far from being as clear as the bulk properties. It seems to be instructive to compare 
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this situation to the boundary behavior of linear phantom chains. The boundary effect for linear random walks can 
be easily taken into account by using the "image" principle which enables to represent the partition function of a 
linear polymer with, say, Dirichlet boundary conditions in terms of a linear combination of two shifted bulk partition 
functions pil |25| | . A priori the same "image" method seems to be inapplicable for randomly branched chains. 

In the present work we generalize the technique developed in the earlier work |l2l for a problem of thermore- 
versible gelation to compute a partition function, Z(N,d), and calculate various thermodynamic characteristics of 
annealed randomly branched loopless N— link polymer in a 3-dimensional space, whose one point is fixed at the dis- 
tance d from the impenetrable wall. The approach used in our work is based on a diagram expansion of a partition 
function of randomly branched polymer in a semi-space. Summing the diagrams we then construct and solve directly 
the Dyson-type equation for the aforementioned partition function. The developed method enables us also to compute 
the density profiles of the monomers having different numbers of neighbours as well as to derive an expression for the 
correlation function, G{ri, r 2 , N), of the branched polymer in a semi-space. 

The paper is organized as follows. In the Section II we present the model of randomly branched polymers and 
discuss in details the methods we are using. The brief review of known main results concerning the randomly 
branched polymers in the bulk is given in the Section III. The Section IV is devoted to the calculation of a generating 
function of branched random walk in a semi-space. That is the central part of our paper. In the Sections V-VII, 
using the expression for the generating function, we compute respectively: the partition function of an iV-link random 
tree, the distribution of branchings in the tree, and the pairwise correlation function. In the last Section we briefly 
summarize and discuss the obtained results. 



II. THE METHODS AND THE MODEL 



In order to make the paper self-consistent, it seems to be instructive to formulate the general thermodynamic 
language which has been elaborated for the description of the thermoreversible gelation process 0> and is very 
convenient for our needs. To begin with, let us consider a system of N identical 3-functional units capable to form 
reversible bonds between each other. The partition function Z{N) of such a system can be written as a product of 
two terms: 

Z{N) = Z str (N)Z int (N) (1) 

where Z stT (N) and Z- mt (N) denote correspondingly the "structural" and "interactional" parts of the partition function. 

The structural contribution to the partition function, i.e. that due to the formation of clusters of specific structure 
reads: 

Z ait {N)= /Et" II [Pafrur^eh-i (2) 

where: 



the product is taken over all pairs of particles forming the manifold with topology T characterized by the 
symmetry index r-r; 

the external sum runs over all possible topological structures T of the system; 
the factor [3 is the weight of pairing; 

the function q(rj, r 1 ) = q(rj — r,) is the probability density to find two connected particles at the points r» and 
Tj correspondingly 0, UH HH • We assume this function to be the normalized Gaussian: 

9(r)=(^Y /2 exv(- S £) (3) 



V27raV 2a V 

where a is a mean-square length of the bond: 

r 2 g(r)dr (4) 



a 2 
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The partition function Zi nt (N) is purely energetic and is due to the interactions among particles and between 
particles and an external field: 

Z int (N)= fe-^^Hdn 

(5) 

U({n}) = Y,^ + zZ v ^ r ^ 

i i,j 

where: 

- the potentials </>( r i) an d Vfri, rj) are correspondingly an external field at the point and a pairwisc interaction 
energy in the system; 

- T is the temperature measured in energetic units. 

In what follows we neglect the volume interactions between particles but keep and therefore we set V(ri,rj) = 0. 

In the present article we restrict ourselves to the investigation of randomly branched polymers only, i.e. we totally 
neglect the possibility of formation of any closed loops of bonds in the system. In that case the calculation of the 
partition function Z(N) becomes rather transparent and can be described by the simple diagrammatic technique (see 
Fig^as an example). 

z(4) = (X X X X) + (X X >— () 

+ (>-( H) + (aA) 

FIG. 1: The diagrammatic representation of the statistical sum Z(N) for N — 4. 



Indeed, let us consider a system consisting of N 3-functional monomers bonded in such a way that no closed loops 
are present. In general, such a system consists of many disconnected trees. It is convenient to describe the structure 
of each tree by a diagram, where each monomer is represented by a vertex (with an assigned radius- vector r), the 
chemical bond between two monomers is represented by a line connecting the vertices corresponding to the bonded 
monomers and if a monomer has less than 3 bonds with other monomers than each missing bond is represented by a 
short line. Next, let us attribute the weights e~^ r ^l T , /3g(n,rj) and 1 to the vertices, long lines (connections), and 
short lines (missing bonds), respectively. Following the prescription @, the weight of each diagram is given by a 
multiplication of all the weights, integration over all space coordinates and, finally, division by the symmetry index r-p. 
Now, the desired partition function Z{N) is a sum running over all (connected and disconnected) weighted diagrams. 
To proceed further let us introduce the generating function 5(A) of all diagrams (connected and disconnected) as 
follows: 

oo 

S(A) = 1 + 

iV=l 

The function 5(A) is just the partition function of the grand canonical ensemble with a fixed chemical potential 
pi = Tin A. As it follows from the so-called first Mayer theorem (see, for example, js3), we can represent 3(A) in the 
following form: 

5(A) = expx(A,#r)) (6) 

where x(A, <j>(r)) is the generating functional of all the connected diagrams (remember that in our approximation all 
connected diagrams are just trees). Accordingly, the partition function of all trees consisting of N monomers is 

Z con (N) = ±-J^dX (7) 
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The simplest way to calculate x(^->4>( r )) is to define first the generating function of all the rooted diagrams (i.e. 
diagrams with one labelled vertex) <r(A,/3, r) = —TSx(\, cj>(r)) / S(f>(r) as it is shown in the equation below: 

T^WWf = \ + X I gir^e-^dn + *£[ fl (r,n) ff (n,n,)n«-* (n)/r *-< 

i=l 



Xe-*(r)/T 6 ' 4 7 * v ' iy 1 4 

/ 6 , ( r ' r i)ff( r > 7 '2)J7e~ 0(ri)/T dr l + ^— / g(r,ri)g(ri,r 2 )g(r 2 ,r 3 ) JJe _0(ri)/T dr-i 

~ / 5(r,r 1 ) ff (r,r 2 ) 3 (r 2 ,r 3 )ne-^)/ T ^ + ^ 1 5 (r, n)fl(r, r 2 ). 9 (r, r 3 ) f[ e -^«)/ r dr. 



(8) 



i=l " i=l 



/ g(r, r\)g(ri, T 2 )g{ri, r%) J^J e~^ Vi ^ T dri + {terms of higher orders in A and /3} 



A^/3_ 
16 



i=i 

The equation (|SJ) can be easily visualized. The corresponding diagrammatic expansion is displayed in Fig[5] 

<r(r, KP)= + X + ~% + ~\ + 

+ V Y +> V - 

FIG. 2: The series of root diagrams for the function er(A,/3, r) up to 4 terms (see Eq.JSJ). 

Note that \~ 1 e^ r ^ T <j{\, (3, r) is a functional of <p(r) and a function of A/3 only, which is due to the fact that 
the number of bonds in a tree- like cluster is always one less then the number of connected units: Nb onc i = N — 1. 
Therefore, in what follows we redefine 

e = A/3 

and rewrite cr(A, /3, r) as 

a{\,e,r)=\e~^/ T Y j C N e N - 1 (9) 

JV=1 

It is noteworthy that the coefficient Ae - ^ 1 *'/ 7 Cn^" 1 has a simple physical meaning — it is equal to the partition 
function of a randomly branched polymer consisting of N + 1 links with one link fixed at the point r. 

Now, it is easy to see that the following relationship holds: 

\d X (\,(3,<t>(r))/d\ = J a(X,r)d 3 r (10) 

Thus, to find the function x( e ) we should just integrate the generating function er(A,/3,r): 

X(A, /3) = ^e"^ £ = W -(A, r, r)dr (11) 

jv=i iV + 1 -Jo 

Since the weights of the different branches of any tree are factorized under assumption of the absence of closed 
loops and volume interactions in the system, the generating function p can be written in the factorized form: 

a(A, £ ,r) = ^ e -^)/^ 3 (e,r). (12) 
o 

Here t(e, r) is a generating function of one branch, which has the following series expansion: 

9 , 2 



t(e,r) = l + ~ J g(r,r 1 )e-^ T dr 1 + e - J g(r, n)g( ri , r 2 ) \{ e-+W T dn 

3 p ^ 3 p ^ 

+ y J g(r,ri)g(r 1 ,r 2 )g(r2,r 3 )l[e'^^ T dn + ~ J g(r, n)g(n, r 2 )g{r u r 3 ) J[ e~^' T dn 
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The diagrammatic form of this expansion is shown in fig. UJl. 



cr(r) 



t(r,e)= = - + — < + 

(a) 

( 

(b) 
(c) 

FIG. 3: Diagrammatic form of: (a) Eq.JTSJ, (b) Eq.JHJ, (c) Eq.fTSl. 

The infinite series l|13|l could be easily evaluated taking into account the aforementioned statistical independence 
of different branches. Indeed, it is easy to see that all the terms in (|13[1 are similar to that of the right-hand side 
of the diagrammatic equation presented in figEJ). Therefore, the generating function t satisfies the following exact 
non-perturbative equation : 

t(e, r) = 1 + | J g(r, n) e -^>/ T t 2 (e, n)dn. (14) 

The equation (|I4|I plays the central role in our paper. Solving it, one gets both "rooted" and "unrooted" generating 
functions p and x an d obtains the partition functions 2 and Z . The method of generating functions enables also 
to compute easily the numbers pf of the monomers having exactly / missing bonds (or ] max — f bonds with other 
monomers). (In what follows we refer to the monomers with / = and / = f ma x — 1 as the junction monomers and 
dead ends, respectively.) The corresponding generation functions are, obviously, 



N=l 



The iV-th term C^ iy) in the series expansion of pf has a clear physical meaning. It is just a partition function 
of a randomly branched (N + I)-link polymer with: i) a root fixed at the given point r, and ii) exactly / branches 
starting from this root. Comparing this definition with that of Cm (see equation |5| and discussion below), we get the 
probability Pfir) for a vertex of a randomly branched ./V + 1-link polymer situated at the point r to be /-functional, 
which equals 

Mr) = cmT (16) 

Now, after the generating function t is calculated, one can also calculate readily the binary correlation function. 
Indeed, the correlation function G{r\ 1 r2) satisfies the following Dyson equation 

G(r 1 ,r 2 )=g(r 1 -r 2 ) + e f g(n - r 3 )t(r 3 ) e -^V T G(r 3 ,r 2 )dr 3 (17) 



which is represented in the diagrammatic form in Fig0] 

Note that equation l|17|l differs from that suggested earlier in |2^| where t(r 3 ) was improperly replaced by 
/ G(r 3 ,r 4 )(ir4. 



III. THE BULK PROPERTIES OF A RANDOMLY BRANCHED POLYMER 



In this section we use the formalism introduced above to re-derive the well-known main characteristics of the 
randomly branched polymers in the infinite space, which will be useful to compare with the results in the semi-space 
to be obtained in the subsequent sections. 



= • — • + •- 



FIG. 4: Visualization of the Dyson equation Eg. (1171 . The correlation function G and the function g are shown by the double 
and ordinary solid lines, respectively 

In the infinite homogenous space (i.e. for <p(r) = 0) the equation defining the generating function of branches t (|14|) 
has a unique solution which is invariant with respect to translations and approaches unity when e — > 0. This solution 
is easy to find as equation (|14|) becomes purely algebraic due to the aforementioned translation invariance: 



t = (18) 
e 

where we took into account that J g(r)dr = 1. Note that all other solutions (in particular that approaching infinity 
when e — > 0) do not have any physical meaning and therefore in case of spatially inhomogeneous external field we will 
be interested only in solutions of ea. ll4|) . which approach the solution given by Eq. I|18f) when r — ► oo; <f>(r) — ► 0. 

The equilibrium densities of junctions and dead ends can be easily calculated in the way prescribed in the previous 
section. Indeed, we have 



,t 3 (2-3e)-(2-e)VT^2e 
P = A 3! =A 3? (19) 

for the generating function of all rooted randomly branched trees. The general expression for randomly branched 
trees with / missing bonds in the root (compare to 1)150) reads 

(t - V\$-f) 

«-»-whv (20) 



Expanding Eqs. 1|19|) . I|2UI) into series with respect of e and substituting the results into eci. (|16[) one obtains finally 
the probabilities of junctions and dead ends: 

(N-2)(N-1) (JV + 5) 

P0{N) = 2(N + l)(2N+iy P2{N) = 2(2JV + 1) (21) 

both values approaching 1/4 as N tends to infinity. 

Now, to calculate the correlation function G{ri,r2) in the bulk we substitute Eq. (|18|l into Eq. I17|l and get 



G(r) = g(r) + (1 - VT^2l) / g(r')G(r - r')d 3 r' (22) 



(note that the correlation function in the bulk depends only on the distance between two roots r — r\ — r%). 

There are various methods of solving l|22|) . we will stick here to the one which seems to be the most suitable to 
be generalized in what follows to the case of the semi-space. As soon as we are interested mostly in the asymptotic 
characteristics of the trees when N is large, we can fully neglect all the short-range (i.e. on the scales of order r ~ a) 
peculiarities of a correlation function. We therefore replace the first term of the r.h.s. of ll'L'l) by a Dirac delta-function 
S(r) (see eqs.©-®) and expand the slowly changing function G(r — r') into the series up to the second order in 
(r — r') (see also |36j). After calculating the integrals we arrive at the following differential equation: 

G(r) = 5(r) + (1 - £) (c(r) + ^AG(r)j (23) 

where we introduced a new variable 



£ = VT^2~f- 
rewrite the 

R = *-^-\/ jzji we can rewrite the equation (|23|l as follows (for R > 0): 



As the function G(r) depends on r only, we can rewrite the Laplace operator as A = r 2 ^ (r 2 -£p)- Denoting 



G"(R) + 2R~ 1 G'(R) - G(R) = 



(24) 
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Therefore, one gets finally the correlation function in the bulk 



G(r) = A e -*±V (25) 



(26) 



which allows us to re-derive the well-known result (see |22j,|23j) for the gyration radius of the tree: 

It is worthwhile also to rewrite the result (|25|l in the cylindric coordinates: 

G(Z, p)=A 6XP { -^ Z2 + p2) = A [ + °° exp (-c\Z\)J (py/J^l) dc 

y/Z 2 + p 2 Ji V / 



(27) 



where Jo(x) is a Bessel function of 0-th order, and the variables (Z, p = ^/ x 2 + y 2 ) are the usual cylindric coordinates 
renormalized by the factor J ■ 



IV. GENERATING FUNCTION t IN A HALF-SPACE 

Assume now that our branched polymer is displaced in a semi-space x > 0. The presence of an impenetrable wall 
situated at x = is described by the potential <fi(r) = <fr(x, y, z), where 

, , * f for x > 

<Kw)H _ fora;<0 (28) 



In this case one can rewrite (|14|l as follows 



c 



t(x,y,z,e) = 1+ - / dx / rfz / dy g(x,y,z;x ,y ,z)t (x ,y ,z ,e) (29) 

To solve Eq. (|29|l we suggest the following procedure. First of all we represent 129fl in the form: 

g- 1 {t-l)= e -f (30) 
where <7 _1 is the inverse operator of <?, the latter being defined as 

gf(r) = J g{r,r')f{r')dr' 
Following [34| we expand (t — 1) in (|3U[1 in the series up to the second order in (x — x'): 

t{x', y', z')-l = t(x, y, z )-l + (x-x')^+ ILL + ... (31) 

where we took into account that for obvious physical reasons t depends neither on y nor on z. This expansion is valid 
for sufficiently smooth function t(r). Therefore, in the case under consideration this substitution is acceptable if we 
are not too close to the surface (the wall). Substituting IpITj) into (J3UJ), we obtain finally 

~ 1 / x , x - a 2 d 2 t(x, e) , „ 

g- 1 (t-l)~t(x,e)-l + — £j (32) 

which results in the following differential equation instead of the integral one: 

/~ n d 2 t(x,e) e 9/ „ . . . 

t(x,e)-l- g ~ 2 ' =-t 2 (x,e) (33) 
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where 



-V6 



is a reduced distance from the surface. 

The differential equation can be solved via the substitution p(t) = M . We get after some algebra: 

dy 



x(e,t) 



to 



3=r+2y2 _ %+ | (l-^-l + 3e 



(34) 



where to = <(0, e) > 1 is a boundary value of t to be specified later (one could not define this boundary condition a 
priori, as the boundary x = does not belong to the region where the substitution (|31[) is valid; see section VI for 
the discussion of the proper choice of t ) and we have already used the boundary conditions at +oo: 



lim t 

x — *+oo 



l - VT^2i 



lim — = 

>+oo dx 

to define the last term in the denominator in l|34fl . 
Performing the substitution 

£ = A /T _ 2i 
a = l + iey-l)^- 1 
we arrive at the simple integral for the function x(e, t): 

l+(et-l)/« 

l+(et -l)/£ 

which we can easily evaluate: 



(35) 



(36) 



(37) 



6| 



~x(U) = 4 farctanl/ 2 + 4 ^ty 2 ^ - arctanhV 2 + * + ~ X > 



Inverting l|38|l we obtain the desired partition function t(x 



t(x) = ■ 

where the auxiliary function <?(£,io) reads 



1 - £ - 3£ sinh" 



V2 + 4e + to(e 2 -i)-yec 



(38) 



(39) 



(40) 



Note, that as x tends to infinity, the function sinh(...) in (|39|l does too and, therefore, the function t(x) approaches 
its bulk value (TH) . 



V. PARTITION FUNCTION OF A RANDOMLY BRANCHED POLYMER NEAR THE SURFACE 



To compute the desired partition function C/v(r) of a single randomly branched polymer, consisting of N + 1 links, 
one of which is fixed at a given point r, we should expand the generating function 



p(X,e,r) 



At 3 (e,r) 



(41) 
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where t is given by i|39|l . in the power series in e. We perform this expansion in two successive steps. First we expand 
t in the power series in £: 



t - 2(1 -a l-3*o 2 + j-^*o 3 5 i + " Id-*) 



63' T ° 3 r° 



5 5 (2 -t ) 5/2 

7/2N 



+ 2 (^) )(l-2 e ) 3 / 2 + 0((l-2 e ) 2 ) 



(42) 



where 



is a reduced coordinate, and we have used l|3tj[l to replace £ by e. Thus, the first singular (with respect to e) term in 
the series expansion of t(e, r) is proportional to (1 — 2e) 3 / 2 and not to (1 — 26) 1 / 2 as in the bulk case. Now, the first 
singular term in the series for p(A, e, r) equals to 

P2 S = ^ *(1 - 3* - 2 f U + 2 (^-) ^ (1 - 2 e ) 3 / 2 (44) 
Now, taking into account that 

3 2 , Q ^ (27V- 5)!! N 



(f-2e) 3 / 2 = l-3e+^ 2 + 3^ 



JV! 

N=3 

and allowing for the asymptotic behavior of the coefficients in the above sum: 

{2N m 5)U = 2 N TV" 5 / 2 (1 + OiN- 1 )) 

we obtain that the desired partition function C n (r) for n» 1 tends to 

C N (x ) = A(x Q )N- e e N (45) 
where the function A(x) depends on the microscopic parameters of the model: 

7/2\ 



A(x ) = | *o 3 (1 - 3* - 2 ) 2 U + 2 (jZ-^) 



(46) 



and the exponent 9 = | is universal. 

If xq is large enough (i.e. we are not too close to the surface), one can neglect all the sub-dominant contributions 
of Xq and can rewrite (|45|) - (|46|l in a simpler form 



N 3 / 2 \Na 4 

Note that the additional (in comparison with the bulk case) factor 



Cn{t) c — - -— (47) 



x 4 



Na 4 yaN 1 /* 

is just the 4-th power of the distance from the surface to the gyration radius of the randomly branched polymer l|26(l . 
Thus, in contrast to the bulk behavior, the partition function of the randomly branched polymer near the surface 
carries the information about the spatial dimension of the branched polymer. 
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VI. THE PROBABILITY OF BRANCHING 



Our next objective is to calculate the probabilities of /-functional branchings, p/(r,N) (/ = 1...3) in a N + 1- 
monomer tree. As we are mostly interested in the case of very large trees we restrict ourselves here to calculation of 
the limiting values Pf(r) = lim pf(r,N). To do that we should find, according to l|16|l . the asymptotic form of the 



N- 



coefficients Cisr(r) and Cjf' '(r) in the series expansions of t 3 (r) for Cjv(r) and of (t(r) — l) 3 , 3(i(r) — l) 2 , 3(i(r) — 1) 
for C$\r), CP{r), C^ 3) (r) respectively. Thus, we find in particular for the probability of junctions (/ = 3) 



C N (r) 



4 




6 








' x o 




•' 


4 




6 




0- 
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)' U + 2 (2^) 2 N N- 5/2 (1 + (iV- 1 )) 



(48) 



and therefore 



P3(r) 



2xq 



(49) 



Similarly, one easily finds for the probability of dead ends (/ = 1): 



Pi(r) — lim 

C N (r) (2x% -&Y 



(50) 



In the Fig. [SJwe have plotted the dependencies pi^(xo). As xq — » 00 these probabilities approach their bulk values 
given by Ea. (|2T)l . 
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15 



20 



FIG. 5: Profiles of junction points, pi, and of dead ends, P3. The dashed line corresponds to 20 = v6. 



On the other hand, if one chooses to, which is the boundary value of the function t(x) to equal unity (1 < t(x) < 2 
for < x < 00 (see ca. (11411 V the value of xq approaches y/6 near the wall and the probabilities of dead ends and 
junctions tend to unity and zero, respectively. This result seems to be rather natural: as the wall is impenetrable at 
x = 0, there should be mostly dead ends. Therefore we arrive to the conclusion that the choice of <o = 1 is at least 
the most natural one thus solving the problem of proper choice of the boundary conditions at x = outlined in the 
discussion after eq. (|34|l . We restrict ourselves to this choice in what follows. 
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VII. THE CORRELATION FUNCTION IN THE SEMI-SPACE 



To find the correlation function G(r 1; r 2 ) of a randomly branched polymer in a semi-space in 3D we have to solve 
the equation (fTTjl where the partition function t(r) given by eq. (|3"§)l . To do that let us multiply fTTf by ei(ri) and 
define a new function 

r(ri,r 2 ) = et(r 1 )G(ri 1 r 2 ). 

After such substitution we get 

r(r 1 ,rs,)=et(ri)( 5 (|r 1 -r 2 |)+ / g (\n - r 3 |) T(r 3 , r 2 )d 3 r 3 ) (51) 



Similarly to what we have done in the previous sections we substitute g{\r\ — r 2 \) = S(n — r 2 ) into the first term of 
the r.h.s. of Q51(l and expand r(r3,r 2 ) up to the second order in (ri — r 3 ). Evaluating the integrals we arrive at the 
following differential equation 

v, p) = e t(u + v) (T(u, v, p) + A u>p T(u, v, p)) (52) 

where 



V&x 2 V&(xi - x 2 ) V 6 ((yi - y?) 2 + ( z i _ z 2) 2 ) ,_„^ 

v = ; u — ; p = — (5d) 

a a a 

In what follows we assume u > 0. Such choice does not lead to any loss of generality due to the symmetry of the 
correlation function: G{r\ 1 r 2 ) = G(r 2 , ri) and the case of u = is to be taken into account via boundary conditions). 

Now we seek for the solution of eq. I|52|) in the form 

T(u,v,p) = J f(c 1 v)T 1 (c,u)T 2 (c,p)dc (54) 

where the arbitrary function f{c,v) is to be determined later by the boundary conditions. The ansatz (|54() allows us 
to separate the variables in (|52|l . We thus obtain a simple equation for T 2 

_! d ( dT 2 \ 2 . . 

1 • ' - c T 2 (55) 



i9p \ dp 
the general solution of which is 

r 2 (c,p) = AJ (cp) + BY (cp) (56) 

where Jo and Yo are the Bessel functions of the first and second kind. Due to the boundary condition r(c, 0) < oo, 
we set B = 0. 

The equation for Ti is as follows: 

<5 2 ri ii , >\„ , , 1 — et(u + v) 

c 2 + <p(u,v))T 1 ; <p(u,v) = \ ' 57 

ou z et(u + v) 

If we plug the exact expression (|3l?jl for t(u + v) into (|F7jl . the resulting equation seems to be not solvable. However, 
if we approximate the genuine function tp{u, v) by its asymptotics in the most interesting regime £ — ► (i.e. N — > oo) 

and a£ -1 / 2 > x\ > a, 

y(u,t;)M , , 2 (58) 

(it + v + pj z 

where p = 2^/6 = lim^^o (£ -1 ^ 2 m g(Q) j the resulting differential equation is solvable (by reduction to the Bessel one) 
precisely. 

The corresponding solution is, however, not too accurate. To obtain an improved solution, which approximates 
the genuine one in the hole range of variables u, v, we are looking for the proper replacement of the genuine function 
f(u, v) defined in l|57|l in the form 
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thus preserving the limiting behavior of tp at u — > oo, £ — > 0. To preserve the limiting behavior of (|57|l at u — > we 
now need to set 

tp(u = 0, v) = tp{u = 0, v) 



dtp(u, v) 



du 



dtp(u, v) 



u=0 



du 



u=0 



which results into 



tp(u) 



I2({p + v) 2 - 12) 
((p + v)(u+p + v) - 12) 2 



(59) 



At last, in the limit of u £ 1/>2 one should replace tp by its bulk value tp = ^/(l — £) which does not depend on u 
and therefore does not affect the desired solubility of the differential equation (|57jl . Thus, finally, 



12((p + v) 2 - 12) 



(60) 



l-£ ' {{p + v)(u + p + v) - 12) 2 
seems to satisfy the desired conditions, both being a good approximation of tp in the whole range of parameters and 

-e t {u-\-v 
t{u+v) 



making the eq. (|57|) solvable. The comparison of the functions tp(u, v) and tp(u, v) = 1 A^ u+ ^ for two different values 



of the parameter v (v = and v = 3) is shown in the FigEI 

1,0- 
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cp 0,1 




FIG. 6: Plots of the functions 
insertions in both plots show the ratio tp/tp. 



tp(u,v) (thin line) and <p(u,v) = ^f^"^^ (bold line) for v = (left) and for v = 3 (right). 



The 



As one sees, the approximation of the exact function tp(u, v) by the function tp(u, v) is reasonable for v — and 
already very good for v = 3. Let us recall that the distance from the wall v = 3 in the non-renormalized (initial) 
coordinates according to (|53|) is X2 = 3a/\/6 ~ 1.22a. 



The resulting differential equation 



r? = (c 2 + m)^ 



is solved precisely via substitution w = u + v + p — ^q^. Taking into account the boundary condition ri — > at 
infinity, one obtains the final result in the following form 

Ti = ^pus K v {c w) (61) 

where K v {dv) is a modified Bessel function of 2-nd kind of order v, and the parameters are as follows: 

y/A9(v+p) 2 - 576 



2{v+p) 



1-C 



(62) 
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Substituting 156(1 , l|61|) into l|54|) we arrive at the desired solution of the equation i|52|) : 

T(u,v,p) = J f(c,v)J (cp)y/pwK 1/ (c'w) dc 
To make this equation comparable with that in the bulk case (|27j) let us introduce the new variables: 



(63) 



C 




(64) 



After such a substitution the correlation function acquires the form 

r(u, v,r) = J f{c, v) mr^c 2 - 1) Vw k v {cw) dc 

The unknown function /(C, V) is to be determined by the boundary conditions 



T(U,V,R) 



tbu\k 



[/»! or V»l 



One thus easily finds 



C[V+P- 
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V + P 



C > 1 



(65) 



(66) 



(67) 



where we have omitted all terms which do not depend on c and v. Let us stress that the condition C > 1 ultimately 
defines the limits of integration in l|67|l . Thus we arrive finally at the following expression for the correlation function 



/+oo 
^CW e c(v + p-i2(v+p)-i) Kuicw) j o ^Ry/cpZTij dC 



(68) 




a) 
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b) 



FIG. 7: a) The correlation functions G(U, V, R) for £ = 1CT 4 , i? = 0.01 and V = 0.05 (dotted line), V = 1 (dashed line) and 
V — 2 (dot-dashed line) as compared to the bulk correlation function for similar R (bold line). All the curves are normalized 
by the condition J G(U, V, R)dU = 1; b) The ratios G(U, V, R)/Ghuik{U, R) for values of £, R and V similar to that in a). 

Figure [7| shows the dependencies G(U) for different values of V as compared with the bulk behavior. One sees 
easily how the presence of the wall affects the correlation function. In particular, it is interesting to note the increase 
of the polymer density in the vicinity of the wall, and also the smaller increase of the density in the outer region due 
to polymer-wall hardcore repulsion. The rapid decrease of the correlation function for V = 0.05 with the growth of 
U is due to the fact that in the system controlled by the fugacity e, the mean size of clusters in the vicinity of the 
wall is much smaller than that in the bulk (this is clearly outlined by the fact that t(x) — > 1 as x — > 0). 
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VIII. DISCUSSION 



Summarizing, in this paper we presented a rigorous procedure to describe the behavior of the large ideal trees near 
the impenetrable (non-adsorbing) wall and obtained rather accurate approximate expressions both for the partition 
function and 2-point correlation functions of the system. 

For comparison of our results with those obtained earlier |32l l3^ | we need a better understanding of the applicability 
of the original approximations of "large" and "no" excluded volume. 

For simplicity, let us fully neglect the attractive interactions between the monomers. In this case the only volume 
interactions in the system are those due to the excluded volume. To estimate how much these interactions disturb 
the original conformation of the non-interacting polymer one should (see [3J) calculate the mean number of pairwise 
contacts of the monomers. This value equals 112 — pN , where p, which is the mean density of the polymer, can be 
calculated as p — Nv/R® ~ (v/a 3 )N 1 / 4 , where v is the excluded volume per monomer. One can estimate that the 
theory presented above should be valid in case when the number of pairwise contacts (and therefore, as one can easily 
see, also the contacts of higher order, i.e. triple, etc.) is negligible, or, in other words, if N -C (a 3 /w) 4 / 5 . On the 
contrary, if N ^> (a 3 /v) 4 ^ 5 , we expect the results of [33| to be valid. Note, that the experimental value of the key 
parameter Li = v/a 3 (we refer to it as the Lifshitz parameter) can vary widely. E.g., for the lattice model Li ~ 1 
and, therefore, practically any tree on a lattice cannot be considered as an ideal one. On the contrary, for many real 
polymer systems Li -C 1. Indeed, imagine, for example, a tree constructed of star-like monomers, whose arms are 
long polymer chains with associating groups at the end. For the polymers of such architecture one can expect the 
Lifshitz parameter Li to be much less then unity. Thus, though in the limit of N — > oo the volume interactions are 
always important, there can be a rather wide range of polymer sizes 1 -C N <C Li~ 4 / 5 where the results obtained 
above are correct. 

As it is mentioned in the Introduction, the statistics of randomly branched ideal polymer chain near the impenetrable 
boundary was studied by supersymmetric methods in the work |33j , where the authors follow the general scheme of the 
supersymmetric " dimensional reduction" for branched random walks formulated for the first time in |26j and exploited 
later in [Hl^i]. The authors 01 the 

paper 33] have computed many thermodynamic properties of branched polymer 
chains near the repulsive and the adsorbing impenetrable surfaces in 3D. For our purposes the most important are 
those of these results, which are related to the "non-adsorbing" (i.e. repulsive) regime. 

In |33| the authors have got the following results below the adsorbing transition point for the quantities of our 
interest: 



G N (z, z') = erfc(C - CO + erfc(C + CO - 2crfc(C) + 16r 2 exp [T(T - 2C - 2C')] 
where 

C = zVn, C = z'Vn, r = eVn 

and E is some constant independent of N. 

It can be seen that the surface critical exponent 6 = | obtained in our work (eqs. Q45|). (|47fl ) coincides with the one 
of |33| (computed also in |32|). We can thus assume that this critical exponent seems to be independent of volume 
interactions, therefore being genuine for all N 3> 1 independently of the Lifshitz parameter Li. 

Now, as far as the correlation functions is concerned, their behaviour depends on the regime significantly. Indeed, 
(|69|l suggests the characteristic length scale of the problem to be of order r ~ y/~N (this result being in compliance 
with the gyration radius for the randomly branched polymers with excluded volume as obtained in 26]). On the 
contrary, our approach results into r ~ N 1 ' 4 , as the integral in Ht)8|) is mostly defined by its lower limit (see also I07J, 
where this characteristic scale is clearly outlined), in full coincidence with that obtained for the characteristic scale 
for the non- interacting randomly branched polymer in the bulk (compare to j2^,|23| and equations l(2*B|) in this 

paper) . 
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